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the elastic axis, equation (56), m 

tensile axis (area centroid) offset from elastic axis, positive 
when in front of the elastic axis, equation (35), m 

Young 1 s modulus , N/ m 2 

shear modulus, N/m 2 

unit vectors associated with undeformed beam coordinate system, 
figures 1-3 
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figures 2, 3 

n fixed n unit vectors with K in direction of rotation of $ in 
jt, figure 1 



blade cross-section moment of inertia from y 1 axis, 
equation (35), m 4 
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equation (56) , m 

blade cross-section integral, equation (56), m 
blade cross-section integral, equation (56), m 

components of generalized nonconservative forces in x, y, and z 
directions, equation (60), N/m 

mass per unit length of the beam, equation (56) , kg/m 

internal structural moment 

twisting moment about the x and x' axes, figure 4, N-m 
S x , - (P x i)'. equation (93), N-m 

bending moment about y and y' axes, figure 4, N-m 
bending moment about z and z' axes, figure 4, N-m 
generalized nonconservative moment about x' axis, N-m/m 
fixed point in £ and £ at the root of the blade, figure 1 
inertial and aerodynamic forces 

warp term in strain energy expression, equation (34) , N-m 2 

inertial and aerodynamic moments 

distance along the deformed elastic axis 

length of blade, m 

inertial frame 

twisting moment arising from shear stress, equation (34), N-m 
time, sec; also blade thickness, m 
tension (same as V x <) , equation (62), N 

kinetic energy, kg-m 2 /sec 
transformation matrix 

a twisting moment arising from longitudinal stress, equation (34), 
N-m 

elastic displacements in the x,y,z directions, respectively, 
figures 2, 3 
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strain energy, N-m 

axial force (tension) in the x, x' directions, respectively, 
N 

shear force in the y, y' directions, respectively, N 
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vector velocity of point P in J? , m/sec; also internal 
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undeformed blade and y toward the leading edge, 
figures 1-3, m 

axis system tangent to the deformed body, figures 2, 3, m 

coordinates of a point in the undeformed blade, m 

coordinates of a point (which was at xq, y Q , z Q in the 
undeformed blade) in the deformed blade, m 


notation for writing the energy expressions in a concise form, 
equations (38) , (53) , and (55) 


precone angle, figure 1, rad 
variation of ( ) 

virtual work of the nonconservative forces, equation (29), 
N-m 
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small parameter, equation (4) 
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NONLINEAR EQUATIONS OF MOTION FOR THE ELASTIC 

BENDING AND TORSION OF TWISTED 

NONUNIFORM ROTOR BLADES 

D. H. Hodges and E. H. Dowell 1 

Ames Research Center 
and 

U. S. Army Air Mobility R§D Laboratory 


SUMMARY 


The equations of motion are developed by two complementary methods, 
Hamilton's principle and the Newtonian method. The former is more mathemati- 
cally rigorous and systematic, while the latter provides more physical insight. 
The resulting equations are valid to second order for long, straight, slender, 
homogeneous, isotropic beams undergoing moderate displacements. The ordering 
scheme is based on the restriction that squares of the bending slopes, the 
torsion deformation, and the chord/radius and thickness/radius ratios are 
negligible with respect to unity. All remaining nonlinear terms are retained. 
The equations are valid for beams with mass centroid axis and area centroid 
(tension) axis offsets from the elastic axis, nonuniform mass and stiffness 
section properties, variable pretwist, and a small precone angle. The strain- 
displacement relations are developed from an exact transformation between the 
deformed and undeformed coordinate systems. These nonlinear relations form an 
important contribution to the final equations. Several nonlinear structural 
and inertial terms in the final equations are identified that can substanti- 
ally influence the aeroelastic stability and response of hingeless helicopter 
rotor blades. 


INTRODUCTION 


The dynamic response and aeroelastic stability of helicopter rotor blades 
can be determined in many cases from linear equations of motion. Much work 
has been done with linear analyses - notably the work of Houbolt and Brooks 
(ref. 1), which contains a systematic derivation of the partial differential 
equations for the coupled bending and torsion of twisted nonuniform blades. 
These equations are not sufficient, however, in cases where nonlinear struc- 
tural or inertial effects are important. Nonlinearities can be important in 
determining the dynamic response of both articulated and hingeless (cantilever 
blades) rotors, and they are especially important in determining the aeroelas- 
tic stability of torsionally flexible hingeless rotors. 


Consultant, Professor, Department of Aerospace and Mechanical Sciences, 
Princeton University. 


This report develops a more complete and general nonlinear theory with 
particular emphasis given to the fundamentals of the nonlinear behavior. The 
theory is intended for application to long, straight, slender, homogeneous, 
isotropic beams with moderate displacements and is accurate to second order 
based on the restriction that squares of bending slopes, twist, t/R, and c/R 
are small with respect to unity. Radial nonuniformities (mass, stiffness, 
twist, etc.), chordwise offsets of the mass centroid and tension axes from the 
elastic axis, precone, and warp of the cross section are included. Other more 
specialized details are not considered, such as blade root feathering flexi- 
bility, torque offset, blade sweep, and droop; nor are configurations 
considered in which the feathering bearing is replaced with a torsional ly 
flexible strap. 

The equations of motion are derived by means of two complementary methods: 
the variational method based on Hamilton's principle, and the Newtonian method 
based on the summation of forces and moments acting on a differential blade 
element. Both methods used together help ensure a more accurate and consist- 
ent treatment of the nonlinear terms. The important nonlinear strain- 
displacement relations, required for both methods, are developed from a clas- 
sical definition of strain and simplified in accordance with the premise of a 
long, slender beam subject to moderate displacements. Applications of the 
equations of motion to rotor blade aeroelastic problems and development of the 
aerodynamic loads are not included in this report. 

Several previous studies have included nonlinearities in the equations of 
motion and identified their importance in forced response and aeroelastic 
stability analyses. A brief discussion of this work should help to acquaint 
the reader with the nonlinearities of interest, the extent to which they are 
treated in previous work, and the rationale for the present development. 

Several approximate analyses have represented the elastic cantilever 
blade with a rigid blade and a spring restrained hinge. Reference 2 showed 
that inertial (Coriolis and centrifugal) nonlinearities were important even if 
only the flap and lead-lag degrees of freedom were retained. Reference 3 (also 
based on the rigid blade representation) showed that structural bending- 
torsion nonlinearities were important when blade torsion was included. The 
flap and lead- lag nonlinear inertial terms were derived for an elastic blade 
in reference 4; however, the radial displacement and an intermediate 
equation for radial force equilibrium were required. 

The extensive treatment of the elastic bending-torsion equations of 
motion in reference 1 retained only linear terms in the final equations, 
although nonlinear terms were included in certain intermediate stages of the 
derivation. In particular, the present equilibrium equations in terms of 
force and moment resultants are identical to those of reference 1. The pres- 
ent inertial terms are also identical to those of reference 1 prior to simpli- 
fication therein. The present work includes a more complete development of 
the strain-displacement relation, and this is required to obtain the elastic 
bending-torsion coupling terms. Reference 5 discusses a simplified form of 
the bending-torsion coupling in the torsion equation - that is, the elastic 
torque due to the product of flap and lead- lag curvatures. The present deri- 
vation includes a more general formulation of these terms as well as the 
corresponding terms in the bending equations. 
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Reference 6 gives a relatively complete derivation (by the Newtonian 
method) of the nonlinear equations for rotor blade bending, torsion and rota- 
tion about flap and lead-lag hinges. Certain second-order nonlinear inertial 
terms due to elastic deformations of the blade are not included, however, and 
these can be important for stability of cantilever blade configurations. 
Reference 7 contains a derivation (also by the Newtonian method) of a system of 
nonlinear equations for the elastic bending and rigid pitching motion of a 
cantilever rotor blade. Most of the nonlinear terms are included, but certain 
inconsistencies in the discarding of higher-order terms resulted in nonself- 
adjoint inertial and structural operators; that is, the stiffness matrix and 
mass matrix are not symmetric, and the gyroscopic matrix is not antisymmetric. 
The nonlinear bending-pitch coupling does not appear explicitly as a struc- 
tural nonlinearity (as in refs. 5, 6, and the present report) since pitching 
moments were derived in terms of the aerodynamic forces. 

Reference 8 presents an earlier and less rigorous version of the 
derivations given in the present report. The earlier derivation was based on 
a simpler development of the strain-displacement relations and several extra- 
neous terms were present in the final equations of motion. Preliminary 
results for the aeroelastic stability of hingeless rotor blades using the 
present equations have been reported in reference 9. 

After a brief discussion of the coordinate systems, transformations, and 
the ordering scheme to be used in this report, the strain-displacement rela- 
tions are developed for an elastic blade. This theory is then simplified to 
second order for long, slender beams. The derivation of the equations of 
motion is divided into two parts. Part I gives the derivation of the equa- 
tions of motion for a rotating beam based on Hamilton’s principle. Part II 
outlines the derivation using the Newtonian method. Both derivations reduce 
to the results of reference 1 when the nonlinear terms are omitted. 

The authors have closely collaborated on the entire report. However, 

Part I is largely the work of the first author and Part II that of the second 
author. Both have benefited substantially from D. A. Peters’ contributions to 
the transformation laws and strain-displacement relations. Numerous useful 
discussions with R. A. Ormiston are also gratefully acknowledged. 


COORDINATE SYSTEMS AND TRANSFORMATIONS 


Several coordinate systems will be used in the present analysis^ The 
orthogonal axes system X, Y, Z and associated unit vectors t, J, K (fig. 1) 
are fixed in an inertial frame Orthogonal axes Xp, y, Z are fixed in a 

reference frame which rotates with respect to $ at constant angular 
velocity ftK. Point 0, a common fixed point of and jB, is located at the 
root of the beam. The plane containing X, xp, Y, and y is called the refer- 
ence plane, or plane of rotation. The x axis, which lies along the elastic 
axis of the undeformed beam, is inclined to the plane of rotation (and to the 
Xp axis) at the precone angle 8p C . The orthogonal axes x, y, z, and the 

corresponding unit vectors T, t, therefore, are also fixed in S. Beam 
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Figure 1.- Undeformed coordinate systems. 

bending deformations shown in figure 2 for 0=0 are described by the 
displacements of the elastic axis u, v, w, parallel to i, j, k, respectively. 
A point on the elastic axis that is located at x, 0, 0 in the x, y, z 
coordinate system before deformation is located at x + u, v, w after 
deformation. 



Figure 2.- Elastic displacements. 










shear center. The cross section is assumed to be symmetric with respect to 
the n axis. The n and X, axes are inclined relative to the y and z axes 
at the "built-in" pitch angle 9(x). When the beam is deformed, the shear 
center for the cross section located at x is displaced an amount u in the 
x direction (not shown) , v in the y direction, and w in the z direc- 
tion. The angle of twist of the cross section changes from 9 about the x 
axis to 0 + <{> about the x’ axis. 

The deformed beam is shown in figure 4 with force and moment resultants 
acting on the face of a cross section. At any point along the deformed beam, 
x' is tangent to the deformed elastic axis. (Note that the distance along 

the deformed elastic axis is also denoted by r, in addition to x'.) The y' 
and z' axes are identical to the n and^C axes, respectively, when the beam 
is deformed. The unit vectors l ' , j ' , k' are parallel to the deformed beam 
coordinate system axes x’, y', z', respectively. Stress resultants and 
moments are subscripted with x', y f , z' to associate them with the deformed 
beam; for example, M x i is a moment about the x' axis. 



Figure 4.- Beam with resultant forces and moments. 

Several coordinate transformations are used, including the simple 

transformation from f, J, | to l, j, A^more difficult transformation 
(appendix) is one from i, J, k to T' , j', k', which is given to second 
order by 

i* = ^1 - i v' 2 - |w' 2 |i + v'j + w'k (la) 
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-+ ■> / I o\ a 

j ' = - [v' cos (9 + <|>) + w' sin(0 + <f>) ] i + I 1 - — v * ]cos (0 + $ + v'w')j 


(l - j w' 2 ^sin(0 + $)k 


(lb) 


k' = - [- v * sin (0 + <fp) + w' cos (0 + <f>)]i - (l - y v ,2 ^sin(0 + $ 


v f w 1 ]) j 


+ 



+ $)k 


(lc) 


where $ (consistent to second order in 


v', w', 4>) is given by 


♦ - ♦ - 



3w 3 2 v 
^ a^ 2 


dx 


The transformation used in equations (1) may be simplified further in many 
(but not all) parts of the analysis by replacing $ by <f> and dropping the 
squares and products of 3v/3x and 3w/3x. The summation of 0 and cj> repre- 
sents the total pitch angle (built-in pitch plus torsion deformation) with 
respect to the x' axis; w is the bending deformation (flapping) in z 
direction; v is the bending deformation (lead-lag) in y direction; and 
3v/3x } 3w/3x are small rotations about the z and y axes, respectively. It 
is convenient to carry 0 + <p through the analysis without making the small - 
angle approximation until the end of the derivation. 


It will be necessary to transform vector components froiri one coordinate 
system to another. For example, a resultant moment vector M may be written 
as 


or 


M = M x i + M J + M z k (2a) 

M = M x » l ' + M y ,j» + M z .k> (2b) 


in terms of components measured in the undeformed or deformed coordinate sys- 
tems, respectively. The relations between M x i, My> , M z » and M x , My, M z are 
determinated by using the unit vector relations of equations (1) and taking 
appropriate dot products of equations (2) . Thus 
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M x = M x i i' • i + My i j ' • i + M z .k' • i 

/ i 2 . 2 \ 

= M x ,(l - ^2” - -~ 2 — J - M y r[v’ cos (0 + <J>) + w' sin(0 + 4>)] 


- M z i[-v' sin(0 + <j>) + w' cos (9 + <J>) ] 


and, similarly, 


where the primed notation is introduced for subsequent use. 


9v 


9w 


V ' = 9x * W 9x 


(3a) 


My = M x .v' + ^1 -^-j[Myi cos(0 + $ + v'w') - M z i sin(0 + $ + v'w') ] (3b) 

M z = M x .w' + ^1 - ^~2~J f M y’ sin ( 0 + ^ + v'w') + M z 1 cos(0 + 4> + v'w')] (3c) 


Other vectors, such as resultant force V and acceleration a, obey the same 
transformation laws. 


ORDERING SCHEME 


In deriving a nonlinear system of equations, it is necessary to neglect 
higher-order terms to avoid overcomplicating the equations of motion. When 
neglecting terms within a large system of equations, care must be exercised to 
ensure that the terms retained constitute self-adjoint structural and inertial 
operators. These self-adjoint operators lead to symmetric stiffness and mass 
matrices and an antisymmetric gyroscopic matrix in the modal equations. 

A systematic self-consistent set of guidelines has been adopted for 
determining which terms to retain and which to ignore. The dimensionless 
axial deflection u/R is generally taken to be of the same order of magnitude 
as the square of v/R or w/R and thus is small with respect to unity. The 
elastic twist ji is a small angle in the sense that sin <)> = <P and cos <J> - 1 • 
The axial coordinate x is of order R and the lateral coordinates are of 
the same order as the chord c, and thickness t, respectively. Both c and t 
are assumed to be at the same order of magnitude as v and w. The warp func- 
tion X is taken to be of order c times t so that the actual warp dis- 
placement will be an order of magnitude less than the axial displacement u. 
These assumptions can be systematized by introducing e, a parameter of order 
v/R or w/R. Hence 


8 


II 

3| 0£ 

0(e 2 ) 

T1 

R 

0(e) 

V 

R ” 

0(e) 

r, _ 

R 

0(e) 

w 

R " 

0(e) 

II 

0 (e 2 ) 

<t> = 

0(e) 

3A/3q 

R 

0(e) 

X 

R “ 

0(1) 

3A/3 ? 
R 

0(e) 


Within the energy expressions, terms of order e 2 are ignored with respect to 
unity. Thus, if the largest terms of the energy expression are OCe 4 ), then 
all terms of 0(e 4 ) are retained (first-order terms), all terms of 0(e 5 ) are 
retained (second-order terms), and generally terms of 0(e 6 ) are discarded. 
There are conditions under which certain 0(e 6 ) terms should be retained. 

These exceptions will be noted below. Since the physical quantities involved 
in the strain energy and the kinetic energy are fundamentally different, a 
scaling parameter for their respective orders of magnitude is introduced in 
the text. This is done so that the ordering scheme may be consistently 
applied to both strain energy and kinetic energy terms. 

In the application of the above guidelines, it is important that the 
ordering be done within the total energy context, or equivalently within the 
virtual work expression (e.g., the w equation times 6w) . Ordering differ- 
ently in one equation than in another without regard to this consideration can 
lead to the above-mentioned symmetry problems, and introduce dissipative and 
circulatory forces into the final equations that are linearized with respect 
to an equilibrium position. Also, the scheme implies that the same order 
terms should be retained in the v, w, <p equations, but terms of one order 
less should be retained in the u equation. 


STRAIN-DISPLACEMENT RELATIONS 


The development of a nonlinear strain-displacement relation is central to 
both the variational and the Newtonian methods for developing a nonlinear sys- 
tem of equations. The use of this relation, together with the generalized 
Hooke's law, permits the strain energy, the force resultants, and the moment 
resultants to be expressed in terms of deformation quantities. Although the 
primary goal of this paper is to develop a second-order nonlinear theory, we 
first develop a general nonlinear strain-displacement relation, which is then 
simplified consistent to second order for deriving the equations of motion. 
This approach is necessary to ensure that a valid approximation for the 
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nonlinear strain-displacement relation is obtained. The general strain- . 
displacement relation would also be useful for future development of a more 
accurate theory. 


Derivation of a General Strain-Displacement Relation 

Several different classical definitions of strain may be found in the 
literature, depending on the mathematical formulation, reference states (based 
on deformed or undeformed positions) , and coordinate systems used. The influ- 
ence of the definition of strain is considered in the present development. In 
particul ar , for nonlinear strain-displacement relations accurate to second 
order, there is no distinction between two commonly used definitions of strain. 

The classic strain tensor e±j in terms of rj and r 0 , the vector 
positions of the same point on the deformed and undeformed blade, respectively , 
may be expressed as (ref. 10) 


drj 



dr 0 



(dr| 

2 [dr dr) dc J [e-^ ] { dn l 

UcJ 


(5) 


where dr, dn, d£ are increments along the deformed elastic axis and two 
cross-sectional axes, respectively. 

The vector position of a generic point on £he^undeformed beam is given by 
(x, y, z) with respect to the unit vectors i, 3, k where (x, 0, 0) is the 
elastic axis. The corresponding point on the deformed beam is given by 
(1) the position of the deformed elastic axis (x + u, v, w) with respect to 
the undeformed axes l, j, k, and (2) the position of the point relative to^ 
the elastic axis [-X(0 + <f>) + , n, c] with respect to deformed axes i', 3*> k'. 
Note in particular that -X (0 + <j>) + , where ( ) + = 3/3r( ), represents an 
axial position where X is the warp function; X(0, 0) = 0. 

Let us now express the above mathematically. After deformation, the 
vector position of a generic point is 


r l 


It 


-t 

1 



( 6 ) 


where the transformation matrix [T] , given in the appendix, relates the 
deformed axes (1', 3', k') and the undeformed axes (i, j» k ) . 

2This derivation is simplified if [T] is retained in abstract form 
throughout. A first-order approximation of [T] introduced at this point could 
result in an incorrect expression for the shear strain component. 
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Before deformation, 


u=v=w=<J>=0 


(7) 


Thus, using equation (6) 



lol + [T] T 


l 

o 

CD 

r 

o 

II 

H* 

<— i. 

5^ 


n 

n 0 



Jo) 

U=V=W=<J)=0 

S oJ 

l_ 


where 


( 8 ) 


^0 


1 U=V=W=<J>=0 


= 


, u=v=w=<f>=0 


X o( T1 0 ’ l u =v=w=<)>=0 


The coordinates n and c, may be regarded as functions of n 0 and 
respectively. 


From equation (A2) , 


[T] 


U=v=w=(j) = 0 


1 0 0 
0 cos 0 sin 0 

0 -sin 0 cos 0 


(9) 


Thus 


r o = U J k J 


n Q cos 0 - s 0 sin 0 
(n 0 sin 0 + £ cos 0 


( 10 ) 


The position vector differentials are given by 


dr 0 = [i 


[ x-ci - M")dr - (x n „ ^ dn * ho ^ d t) 9 ' 


1 , dn« d£« 

-0'x (n 0 sin 0 + c °s 6)dr + cos 0 dn- sin 0 dt; 

dnn d^ 0 

0'x + (n o cos 0-^o s ^ n dr + sin 0 dn + -^r cos 0 dc 


( 11 ) 
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'.'1 

dr 1 = [ i * 

j ' 

k '] 

[ T ]< 

i 


■dr + [T] [T] 


T + 


|-X(0 + ♦) 
n 



f-XC© + <f>) ++ dr - (X n dn + d 0 (0 + <!>) +, | 


+ 4 

d " 

► 


( d ? J 

— 


( 12 ) 


where 


_ <bc 
" dr ; 


. _ 3X . > - £1 . 

A n " 3n 5 C ' 9c ’ 




3C 0 


Note that 


[T][T] T = 



Thus 


[T][T] T+ + 


[T] + [T] T = 



and 


[T][T] T+ = -[T] + [T] T = 


-U) k 

0 


u -OJj 


(Hi 


w-: 


- 0 ). 

l 


(13) 


where explicit expressions for u)i, wj , to k are given in the appendix. Hence 
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(x + + u + 1 


“ 0 


“j" 


f-x(e + <t>) + \ 

dr : = [V 

y 

k'j 

[T ]•< 

v + 

> dr + 

“k 

0 

"“i 

i 

- 




- 

L W + i 


-0) • 

L J 

u>i 

0 _ 


? I 


►dr 


^-A(0 + 4>) dr - (A n dn + A^ d£) (0 + <|)) + 
+ < dri 

dC 


( 14 ) 


The expressions for the position vector derivatives can be simplified by- 
eliminating x + in equation (11) and [T] in equation (14). On the elastic 
axis n = dn = C = dc = A = A n = A^ •= 0. Thus, 


'x + + u +1 


dr. 


dr 


= Li' r k 1 J [T]< 


n=c=o 


w 


(15) 


However, dr : /dr at n = £ = 0 must be i’ by definition of the i', j', k' 
system. Therefore , 


Li- j- 


k'j 


il 

I 

( x + + u + \ 



1 

0 

i — i 

•i— > 

*H 

II 

V+ 


1 

l „♦ j 


(16) 


The first component of this vector equation yields: 


/l - v +2 - W +2 (x + + U + ) + V +2 + W +2 = 1 


Hence, 


= /l - v +2 - w +2 - u + 


(17) 
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Equation (14) for dry may be now simplified, by using equation (16), to 
yield 






'dr) 


0 

-0) k 

“j" 


r-x(o + $) + ' 








0 



1 

drj = Li' 

V 

k'j 

i 

0 

► + 

“k 

-“i 

< 

n 





,0 , 


_-0)j 

U)i 

0 _ 


« J 


f-A(0 + <)>) ++ dr - (X n dn + d?) (9 + <J>) 

dn 
dC 


+ v-i 


(18) 


With equation (17), equation (11) for 


dr 0 becomes 


dr n = 


U 3 


> 

(/l - v +2 - w+ 2 - u + ) (1 - A 0 6")dr - 9* (^r\ 0 dr) + X Z 0 d? ) 

-6<(/l - v + 2 - W +2 - u + ) (n 0 sin 9 + C 0 cos 0)dr + ^- cos 0 dn - -^-sin 0 d^ 




) UT1 q d£ q 

(n 0 cos 0 -C o sin 0 )dr + -^sin 0 dn + ^cos 0 d? 

(19) 

Equations (18) and (19) can now be substituted into equation (5) to solve 
for the components of the strain tensor cyj . 


2e n = (-n»k + ^j) 2 + |> k (0 + <f>) + + ^il 2 + [Auj (0 + <J>) + + nui] : 


+X 2 (0 + 4>) ++2 + 2(-nu k + Cui-j) - 2X (0 + <|>) ++ (1 - nw k + ?mj) 

+ v + 2 + w +2 - u +2 + 2u + /l - v +2 - w +2 - x +2 0 • 2 (n Q 2 + C 0 2 ) 

+x +2 (2X o 0" - X o 2 0" 2 ) ( 20a ) 
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(20b) 


2ej2 ^nCS + *10 (9 + 40 + “ ^^(9 + <fO + (l - tiojj^ + Cojj) - [Ao) k (0 + 4>) + + £i0jj 
+x+0 ’ ^noa-^e’O + 5 0 ] 


2ei3 - AA^(0 + <f>) (0 +■ <f> ) + - A^(0 + <J>) + (1 - na)^ + £u) j ) + [Au)j (0 + <j>) + + nw^] 


+x+6 ' ds° [ x ?o^-M")-n 0 ] 

(20c) 

2 =22 - 1 -( dn °) * yce* «* 2 - »n„ e ’ 2 (inf) 

(20d) 

2c 23 = . «* 2 - X no *, o (^) S' 2 

(20e) 

2e 33 - 1 - (^) ♦ - x 2 0 6' 2 

(20f) 


where 

0 + = 0'X + 0 ++ = 0"x +2 + 0'x ++ 


Equations (20a) - (20£) are the classical strain components in terms of 
displacements u, v, w, and <f>; the warp function A; the curvatures cu^, w-:, 
and u k ; the cross-section coordinates n, n 0 > S, and £ ; and the pretwist J 0. 
The warp function is determined through Laplace's equation for the cross sec- 
tion. The curvatures and twist are known (appendix). A relationship between 
n 0 , C 0 > an< l o> C is needed to express the strain components entirely in terms 
of displacement variables. This relationship will only be determined for the 
second-order approximate theory. 


Second-Order Approximation for the Strain Components 

Equations (20) form the basis for a general nonlinear strain-displacement 
theory. This theory will now be reduced to second order. First, dn 0 /dn and 
d?o/d£ must be expressed in terms of the elastic displacements and the known 
parameters. The assumption of uniaxial stress, valid for long slender beams, 
( ct 22 = °33 = & 2 3 = 0) and Hooke's law gives 


e 2 2 = -venj 


e 3 3 = -venj 

(21) 


15 


Equations (20d) and (20f) are now 


- 2 » e „ = i - * »> * «* 2 

■2vei i = 1 - + ( e + *) + 


and solving for (dri 0 / dn3 and 



- 4 0 o 

(d? 0 /d?) gives 

1 + 2ve n + X n 2 (0 + <J>) +2 
2 2 
i + S 9 

= 1 (to second order) 

9 

1 + 2ve 11 + X^ 2 (0 + 4>) + 

_ - 2 ” 

1 + X t 0' 

= 1 (to second order) 



( 22 ) 


(23a) 


(23b) 


Thus to second order n and £ are equivalent to n Q and £ 0 , respectively, and 
X = X 0 . 

The strain components z n , z 12 , and c 13 in equations (20) may now be 
reduced to second order. The expressions for (dn 0 /dn) and (d5 0 /dC) are sub- 
stituted in equations (20), higher-order terms are eliminated and ( ) is 
substituted for ( ) + since dx/dr = 1 to second order (eq. (17)). 


: n = u' +^ + ^-Xf'+ (n 2 + C 2 )^'*' + 

_v"[n cos(9 + <(>)-£ sin(0 + 40] - w"[n sin(0 + 40 + S cos(0 + 4>) ] C 24 ) 


e 12 


- \ (c + V*' 


(25) 


e 13 = 2 (o ” 


(26) 


Although the 4>' 2 term in equation (24) is formally negligible, since it is 
of 0(e 4 ) compared to other terms of 0(e 2 ), it leads to a tension term pre- 
viously identified in reference 1 that contributes to the elastic torque T x < 
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about the x' axis (see eqs. ( 33 ) and ( 34 ) in Part I). Hence, it is retained 
here although for certain applications it could be dropped without loss of 
accuracy. 

The shear strain component e 23 is now considered. Equation ( 20 e) shows 
that e 23 is not identically zero, contradicting the previous assumption that 
o 23 7 0. However, e 23 is two orders of magnitude smaller than the shear 
strains e 12 and e 13 . Therefore, neglecting e 23 and hence o 23 is well 
justified. This completes the reduction to second order of the classical 
strain components in equations (20). 

The relationship between the simplified classical strain tensor and 
engineering strain is now discussed. Since the uniaxial stress assumption has 
been invoked (o 22 = a 23 = a 33 = 0), the expression for strain energy in terms 
of the classical strain components is equivalent to the more familiar form of 
strain energy based on engineering strain (the latter form is used in Part I). 
Thus, the engineering strain components can be written in terms of the 
simplified classical strain components. 


e xx = 
e xn = 
e x£ 



( 27 ) 


There are other possible definitions of strain that could have been used 
in this development. Here the strain has been defined on the basis of an 
increment of deformed length dr. The most frequently used alternative 
employs the original length dx. By appropriate transformation, the two 
definitions can be compared. For example, examining the case of longitudinal 

displacement only, u/0, v = w = <j> = 0, the present definition of strain 
gives 


e ll 


du _ 1_ /du\ 2 
dr 2 \dr/ 


From the definition of r 


( 28 a) 


and 


dr 

dx 


1 


+ 


du 

dx 


du 

dr 


1 


du 

dx 
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Thus, the alternative definition of strain gives 



which is the more familiar result. To second order, however, the two results 
are equivalent. 
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PART I 


DERIVATION OF EQUATIONS BY HAMILTON'S PRINCIPLE 


The equations of motion and boundary conditions for a cantilever beam 
rotating at constant speed are obtained from Hamilton's principle. These 
equations are valid for rotor blades that can be represented by a long, 
slender, homogeneous, isotropic beam. This representation includes spanwise 
variations in the mass and stiffness properties, variable mass centroid axis 
and area centroid axis offsets from the elastic axis, variable built-in axial 
twist, and a small precone angle. The external forces acting on the beam are 
characterized by a set of generalized distributed loads. The use of an actual 
helicopter rotor blade as a model would be a more formidable task because of 
the nonhomogeneous , anisotropic structures found in typical blades. The 
present structural representation is similar to that of reference 1. 

Hamilton's principle may be expressed as 

t 2 

[S(U - T) - 6W]dt = 0 (29) 

1 

where U is the strain energy, T is the kinetic energy, and SW is the 
virtual work of the external forces. Suitable expressions for SU, ST, and 6W 
are now determined and then combined to give the desired equations. 


X 


Strain Energy Contributions 

The usual expression for strain energy in terms of engineering stresses 
and strains, defined in equation (27) in terms of the classical strains, is 



+ a x£ e xs)^ n 


(30) 


The first variation is 


&U 


X R XX 

* o •'•'a 


( a xx^ e xx + 


°xn <5e xn 


a X f;6ex^)dn d£ dx 


(31) 


where 


a xx 

^ e xx 

a„„ 

= G G 

xn 

xn 

a x£ 

= Ge xs 
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and 


< 5 e xx = (5u' + v'6v' + w'Sw' + (n 2 + C 2 ) (0 + <)>) ' ' - X5(f>" 

-[n cos(e + <JO - Z sin(0 + <J>)] (6v" + w"6<|)) 

-[n sin(6 + <)0 + ? cos(0 + ij>) ] (6w" - v"6<t>) 

6e xn = -C 6 4> ' 

6e x ^ = n6 4> T - 


( 32 ) 


where we introduce the notation n=n-^^;C = ? + ^rr Since the strains ^ 
(and hence the stresses) are composed of the sum of terms of order e and e , 
any product of stress and strain consists of terms of 0(e 4 ), 0(e b ), and 0(e ). 
Consistent with the ordering assumption (e 2 << 1), 0(e ) terms are neglected 
because of the presence of 0(e 4 ) terms. In terms of stress resultants and 
moments, the strain energy variation becomes 

A R 

6 u = I { V x t (6u* + v'6v' + w'6w') + (S x i + T x » ) 6<f) ' + P x i6<f>" 

+ [M z t cos(0 + 4>) + Myi sin(0 + <|>)](6v" + w"6<|>) 

+ [M z i sin (0 + 0) - M y i cos (0 + <|>)](<5w" - v"6<)>))dx (33) 


where the stress resultants and moments are defined by 

Vx , = JJo xx dn dt; = Ea|u' +^- + ^+k A 2 0-(f.' -e A [v" cos(0 + <fr) + w" sin(6 + <f>)]| 

V = /J>x, - ^ a xn^ dn dc = GJ *' 

MM /v T ^W ,2 \ 

T X ' = JJ (0 + <J>) 'o xx (n 2 + S 2 )dn dz « EAk A 2 (6 + 4) ' (u. + _ 2 _ + ~~ 2 ~ J 

+ EB^e 12 ^ 1 - EB 2 *0' (v” cos 0+w" sin 0) 
p x . = -JJ>xx dn dz = EC^" + EC 1 *(w" cos 0 - v" sin 0) 

My, = // ro X x dn dz = EIy.[v" sin (0 + 4>) - w" cos(0 + ♦)] - ECx*^" 


M z • 


= -JJ na xx dn d£ = EI Z ' [v" cos (0 + *) + w" sin(0 + ♦) ] 


/ , V- 2 w' 2 \ 

EAe A (u’ + — + —J - 


eb 2 *6 ' 4> 1 


(34) 
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Three third-order linear terms (doubly underlined) are retained in these 
equations although, consistent with the present ordering scheme, they could be 
neglected. These terms (and two subsequent inertial terms) contribute to the 
final torsion equation. In the special case of rotor blade configurations of 
very low torsional rigidity, these terms may contribute substantially to the 
magnitude of the torsion natural frequency; hence, they will be retained. The 
linear term ECi<)>" in P x t is usually neglected for beams of closed cross 
section. 

The section integrals in equations (34) are defined as follows: 


c* 

«-» 

III 

< 

dZ 

Ae A 

= J 7 A 1 dn dc 

HH 

X 

III 

vTi 

N> 

dn dz 

I Z ' 

= ffpjl 2 dn dz 

CM 
CT 
v — ' 

III 

CM 

< 

r* 

< 

+ Z 2 )dr) dz 

J 

= SS A (r \ 2 + ; 2 )dn dc 

Bi* = J/ A (n 2 

+ z 2 ) 2 dn dz 

B 2 * 

= // A n(n 2 + c 2 )dn dz 

Ci = // a A 2 dn dz 

Cl* 

- // A ? A dn dz 


These integrals are to be evaluated only over the portion of the blade cross 
section that is structurally effective. The blade cross-section area effec- 
tive in carrying tension is A; Iy> and I z t are the flapwise and chordwise 
moments of inertia, respectively; Ak A 2 is the polar moment of inertia; and J 
is the torsional constant including cross-section warping. Without warping, 

J is equal to the polar moment of inertia. The tension axis offset from the 
elastic axis is e A . is the warping rigidity, and Bj* and B 2 * are sec- 

tion constants equivalent to those found in reference 1. Cj* is not included 
in previous references. 

Because of the antisymmetric character of X and the assumed symmetry of 
the cross section about the n axis, the following integrals involving n, Z 
are set equal to zero and hence do not appear in equation (34) : 


J7 a A dn d? = 0 

ff A Z dn dz = 0 

SI A Z( n 2 + C 2 )dn dc = 0 


/J^A(n 2 + C 2 )dn dz = 0 ' 
J7 a nC dn dz = 0 > 
dn dz = 0 J 


(36) 


These integrals would be nonzero if the cross section were asymmetric. 
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Integration by parts of the strain energy yields 


&U = / 0 R (Y u 6u + Y y 6v + Y w 6w + Y^5*)dx + b(tf) 


( 37 ) 


where 


Y u = “ (V x ' ) ' 


Y v = [M z , cos(6 + <|>)+ My. sin(0 + - (V x .v’)' 


Y w = [M z . sin(0 + <{>)- My. cos(0 + <)>)]" - (V x ,w ')' 

= (P x ,)" - (S x - + T x 0 ' ~ v"[M z . sin(0 + <f>) - M y . cos(0 + <t0 ] 
-w" [M z » cos(0 + <f>) + M y t sin(0 + <f>)] 


(38) 


and 


R 


b(tO =V x .6u|„ + {V x »v' - [M z . cos(0 + <t>) + M y > sin(0 + <|>) ] 1 )<$v | 0 


R 


+ [M z i cos (6 + ij>) + My i sin(0 + <|))]6v' |q + t M z' sin(0 + <(>)- M y . cos(0 + <J>) ]6w' | Q 


+ {V x iw' - [M . sin(0 + <J>) - M y . cos(0 + 4>)] ' )<Sw|^+ [S x . + T x . - (P x O ']H 


iR 


+Pv i 6 <}> 


X TFT L Ai Z 

. I R 


(39) 


It is generally essential to include b (U) in equation (37) when using 
modal solution methods, as discussed below. 


Kinetic Energy Contributions 

The position of an arbitrary point after the beam has deformed is given 
by (xi, y 15 zi) where 

xi = x + u - X<t>' - v' [n cos (0 + <f>) - C sin (0 + <(•)]- w' [n sin (8 + <)>) + £ cos(0 + <f>) ] 

yi =v + n cos ( 0 + <)>)-? sin(0 + <(>) 

zi = w + n sin(0 + $) + C cos(0 + 40 
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or 


X! = x + u - A*» - v'Cyj - v) - w 1 ( z j - w) 
y i = v + Cyj - v) 

Zj = W + (zi - w) 


. (41) 


where -A<f>' is an axial displacement due to warp. The velocity of this point 
on the blade with respect to ./?, the inertial frame, is 


V 



RK x r 


(42) 


, -f 

where r ^Xji^ y } j + zik, and 6/6 1 is the derivative in the rotating frame 
M. The fiK x r term is the velocity contributed by the rotating coordinate 
system. The velocity in $ is 


6r 

6t 


xji + y 1 j + 


(43) 


and the rotating coordinate system contribution is 


fiK x r - -Qy 2 cos B pc i + (Rxi cos B pc - Rz! sin $ pc ) J + fiy 2 sin 6 pc l< (44) 
Thus, the total velocity V is given by 

V = (Xi - S2y 2 cos g pc )i + (y t + fixj cos B pc - ftzj sin 6 pc ) j + (zi + Ry 2 sin B pc )£ (45) 
The kinetic energy T is the volume integral of (l/2)pV-V, or 

T = \ / 0 R /Sa^ • ^ d n d? dx (46) 

and the variation is simply 

ST = r R // A p V • 6 ^ dn dc dx (47) 

where 

6^ = (6xj - R cos Bpcdyj)! + (6y 2 + R cos B pc 6xi 

-R sin B pc 6zi)j + (6zx + R sin B pc 6yj)^ (48) 
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The integrand of equation (46) then becomes 

% • 6^=x 1 6x 1 -nx 1 cos - fiyi cos B pc 6xx + fl 2 yj cos B pc 6yx 

• +Q 2 x ! cos 2 3 pc 5xi - n 2 xi sin B pc cos B pc 6zx + fbcx cos Bp C 6yx 

-n 2 zx sin Bp C cos Bp C <5xx + fi 2 zj sin 2 Bp C 6zi - fizi sin Bpc^x + cos Bp C <$Xx 
-fiyx sin Bp C 6zi + yi<5yj+n 2 y-^ sin 2 BpcSy^ + ^yj sin Bp C <Szx 

+QZi sin B pc <5y 1 + zx<Szx (49) 

According to the variational method, equation (49) must be integrated in time 
between two arbitrary points in time, tx and t 2 • The initial and final values 

(e.g., iqSXil^ 2 ) are taken as zero. Hence, we may anticipate integrating by 
parts and combine various terms in equations (47) and (49) to obtain 

6T = / 0 R /J^[(-Xx + 2Ry x cos 6 pc + fi 2 x 1 cos 2 B pc -^ 2 zi sin B pc cos B pc ) < 5x 1 
+ (n 2 yi - yj - 2f2xx cos Bpc + 2Qzx sin Bpc^^yi 

+ (-ft 2 xi sin B pc cos B pc + n 2 zj sin 2 B pc - 2ftyx sin B pc - z x ) 6 zx]p dn dc dx 

(50) 


where 

Xx = u - (v 1 + 4>W) (y T _v ) - (”' - *V') (Zx -w) - X<f> ' 
yx = v - ( z x - w) 


zx = w + Kyi - v ) 

Xx S u - (v' + $w' + 24>w') (yx - v) - (w' - $v' - 24 >v') (zx - w) - XK 


yx = V - <Kzx - w ) 
zx = w + Kyx - v) 


(51) 


Sxi = 5u - (yx - v) (6v ' + w'6<)>) - (z 2 - w) (6w' - v'5<(.) - X6<t>' 


6yx = 6v - (zx - w)6<t> 


6zx = 6w + (yx - v) 6<f> 
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These expressions have been truncated consistent with the ordering scheme 
(e 2 << 1). Assuming that 6p C = 0(e) the variation of the kinetic energy 
becomes 


ST = / 0 R *^a (Zu6u + Zy6v + Zw6w + z 


+Z v i6v' + Z wf 6w' + Z^t 6<|) , )p dn d£ dx 


( 52 ) 


where 


Z u = fl 2 x + 2Q.v 

Z v = G 2 [v + (y 2 - v)] - v + $(zj - w) + 2^6 pc w 

-2fi[u - v'Cyj - v) - w’Oj - w)] 

Z w = -ft 2 B pc x - 2fi6 pc v - w - - v) 

Z^ = -fi 2 xfyi - v) w 1 + fi 2 x(z 1 - w) v ' - fl 2 [v + (y : - v)](z! - w) 

+v(zj - w) - w(y 2 - v) - n 2 gp C xCy 1 - v) - $[(yi - v) 2 + ( Zj 
Zy» = - (ft 2 x + 2fiv) Cy i - v) 

Z^i = - (ft 2 x + 2ftv) (Z} - w) 

Z^r = -Q 2 Xx 


} (53) 


w) 2 ] 
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The doubly underlined higher-order term could be dropped consistent with the 
ordering scheme, but this would eliminate torsion inertia from the torsion 
equation. As discussed previously, this term is important for small values of 
torsion rigidity and hence is retained. 

Integrating over the blade cross section, equation (52) becomes 
(ST = J (Z u 6u + Zy6v + Z w 6w + Z^Sif) 

+Z v ,Sv f + Z w .6w')dx (54) 
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where 


Z u = m(fi 2 x + 2ftv) 

Z v = m£2 2 [v + e cos(0 + 40] + 2mQ(8p C w - u) 

+2meft(v' cos 6 + w ' sin 0) - mv + me^ sin 0 
Z w = -m8 pc (ft 2 x + 2fiv) - mw - me$ cos 0 
Z* = -mk m 2 $ - mfi 2 (km 2 - km^cos (0 + f)sin(0 + f) 
-mefi 2 x(w' cos 0 - v' sin 0 ) - meft 2 v sin 0 
-meQ z 8 pc x cos 0 + me(v sin 0 - vi cos 0) 

Zv* = -me [fi 2 x cos(0 + 40 + 2Uv cos 9] 

Z w , = -me [Jl 2 x sin(9 + 40 + 2-Ox sin 0] 


(55) 


and where the sectional integrals are defined by 

ie // A p dn dc ; me=/J^pn d n d£ 

mk^ £ // A P? 2 dn d? ; mk^ 2 ee J*^ P n 2 dn dC > (5- 

k^ + ^m 2 = k m 2 ; d n d C = dn = JX pA dn d ^ = 0 > 

assuming cross-section symmetry about the n axis and an antisymmetric warp 
function A. The terms involving Cy 1 - v ) and ( z\ - w) , introduced for 
convenience in equation (41) , are given by 

/Ja p ( yi " V ^ dn = me cos ^ 0 + ^ ' 

// A p( Zl - w)dn dc = me sin(0 + 40 

► (5 

ffpp ( y i - v)( Zl - W)dn dc = m(km 2 - k^Jcosfe + <j>) sin (9 + 40 

JXp[(yi v ) 2 + ( Z 1 “ w) 2 ]dn dc = mk m 2 

A. / 
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The blade mass per unit length is m; the center of mass offset from the 
elastic axis is e; the polar mass moment of inertia is mk m 2 . 

After integration by parts, the variation of the kinetic energy becomes 
ST =J^ R [Z u <5u + (Z v - ZvOSv + (Z w - ZwO«w + Z^Jdx + b(:T) (58) 


where 


b(T) = Z v ,6v| R + Z w .<5w| R (59) 

The kinetic energy contributions are equivalent to the inertial terms of 
Part II. The inertial terms of reference 1, before simplifications were made 
therein (eqs. B. 10-15 of ref. 1), are also equivalent to the present results 
for Bp C = 0. 


Generalized Nonconservative Forces 

The virtual work 6W of the nonconservative forces may be expressed as 


6W 



(L u <5u + L v 6v + L w 6w + M^Stf^dx 


(60) 


where L u , L v , L w , and are the distributed loads that act in the x, y, 
and z directions and a twisting moment about the elastic axis, respectively. 
The nonconservative forces include all forces other than the elastic and 
inertial forces (e.g., aerodynamic or mechanical). 


Summary of Partial Differential Equations of Motion 
and Boundary Conditions 

According to Hamilton’s principle (eq. (29)), the strain energy and 
kinetic energy must be combined into one variational statement. Thus, the 
relative order of magnitude of the inertial terms with respect to the struc- 
tural terms must be considered. The ratio of inertial terms to structural 
terms can be expressed as a dimensionless scaling parameter whose order of 
magnitude is assumed to be 


mfi 2 R 2 

EA 


0(e 2 ) 


Hence, the order of terms derived from the kinetic energy is, in effect, 
increased by 0(c 2 ) so that the largest terms of the kinetic energy that are 
0(e 2 ) become ()(£**), the same order of magnitude as the largest strain energy 
terms. Note that the above scaling parameter is essentially the square of the 
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ratio of the tip speed to the speed of sound in the blade material. For other 
applications, the order of magnitude of this scaling parameter may be modified 
accordingly. Now, for the total energy expression, the first-order terms in 
either the kinetic energy or strain energy of 0(e) are retained along with 
the second-order terms of 0(e 5 ). In general, the third-order terms of 0(e) 
are discarded except for certain third-order linear terms as noted above. The 
variational operations can now be performed on the energy expressions as they 
are written with no further terms discarded. 

By substituting equations (37) , (58) , and (60) into equation (29) and 
using equations (34), (38), and (55), one can obtain the total variational 
equation in terms of u, v, w, and For arbitrary, admissible variations 
6u, 5v, 6w, 6<(>, the coefficients of the variations must vanish in the inte- 
grand for all x from 0 to R and also must vanish in the remaining terms 
evaluated at 0 and R. The former condition will yield four nonlinear partial 
differential equations for u, v, w, and <f>. The latter will specify the 
boundary conditions at the ends of the beam. The four equations are as 
follows: 

<5u equation: 

-T* - m(fi 2 x + 2flv) = L u C61a) 


6v equation: 

-(TV)' + |-EAe A ^u' )c° s (e + *)- EB 2 *0'<f>' cos 6 - ECi**" sin 6 

+ [EI Z * cos 2 (0 + <(>)+ Ely* sin 2 (0 + 40 ]v" + (EI Z * - Ely* ) cos (0 + j>)sin(0 + 4>)w"j 
+ 2nftu + mv - me$ sin 0 - 2m e0(v' cos 0+w' sin 0)-mfl 2 [v + e cos(0 + 4)] 

-2mfi3p c w - {me[n 2 x cos(0 + <t>) + 2^v cos 0]} ! = L v (61b) 


6w equation: 


-(Twl)' + |-EAe A (u’ + sin(e + ♦) -EB 2 *0'<f>' sin 0 + ECj**" cos 0 

+ (EI z f -EI y .)cos(0 + 4 )sin( 0 +i)v” + [EI z . sin 2 (0 + *) + EI y * cos 2 (6 + ♦) ]w”J 

+ mw + mecji cos 0 + 2mftBp C v -{me[Sl 2 x sin(0 + <J>) + 2Q.V sin 0] } ' = L„ - mfi 2 8p C x (61c) 
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5 4> equation: 


EAk A 2 (0 + <f0 



+ EB!*0 ' 2 <j> • 


n » 


EB 2 * 0 ' (v" cos 9 + w" sin 0) 


-EAe A 




(w" cos 0 - v" sin 0) - (GJ<J>') ' + [EC^" + EC 1 *(w" cos 0-v"sin 0)]" 


+ CEI Z » - ElyO [ (w" 2 -v" 2 )cos 0 sin 0 + v"w" cos 20] + mkn^+ mft 2 <|>(km 2 - k^^cos 20 


+me[^ 2 x(w l cos 0 - v' sin 0) - (v - ft 2 v)sin 0 + w cos 6] 


2 2 

M<j> - mfi 2 (k m? - k mi ) cos 0 sin 0 


mi- 


men 2 3p C x cos 0 


(61d) 


where 


V x . = EA 


u' 


w i 2 

+ k/^ 2 0 ’ <]> ' - e A [v" cos(6 + <J>) + w" sin(0 + <|>)] 


(62) 


The boundary conditions, from equations (29), (37), and (58), become 


b(i') - b(T) = 0 


(63) 


Equations (34), (39), (55), and (59) may be used to express the boundary 
conditions in terms of the deformation: 


29 



b(£/) - b(T) = T 6 u| q + ^Tv' - |[EI z i cos 2 (0 + £) + EI y . sin 2 (0 + 0) ] v" 

/ v ' 2 w' 2 \ 

+ (EI z t - Ely.) cos (0 + £)sin(0 + 0)w" - EAe A fu' + ^- + Jcos (0 + 0) 
-EB 2 *0'0' cos 0 - EC^ *0" sin ©j + me[Q 2 x cos(0 + <(>) + 2ftv cos 0]^ 6 v| g 
+ j[EI z . cos 2 (0 + 0) + Ely i sin 2 (0 + 0) ]v" 

/ v’ 2 w ,2 \ 

+ (El z . - Ely.)cos(0 + f)sin(0 + 0)w" - EAe A lu' + ~j~+ -y ) cos (0 + 0) 
-EB 2 *9 1 0 ' cos 0 - ECj*0" sin ojfiv' 

+ ^Tw ' - |(E I z i - Ely. )cos (0 + 0) sin (0 + 0)v" + [EI Z . sin 2 (0 + 0) 

+EIy. cos 2 (0 + 0 ^)]w" - EAe A ^u' +^ 2 ~ +^ 2 ~ ^sin(0 + 0^)- EB 2 *0'0' sin 0 
+EC 1 *0" cos 0 J + me[fi 2 x sin (0 + 0 )+ 2fiv sin 0]^<5w|^ 

+ j(EI z . - Ely*) cos ( 0 + £)sin(0 + 0)v" + [EI Z * sin 2 (0 + 0) 

/ V 1 2 W 1 2 \ 

+EIy i cos 2 (0 + j>) ]w n - EAe^fu T + *^~ + + 40 

3 } 6w ' i o 


-EB 2 *0'0' sin 0 + EC^" cos 0 


+ IGJ0' + EAk A 2 (0 + 


£) ' 


v» 2 W ’ 2 \ 2 

' +X 2“ +} V) + EB i* e ' ♦ 


-EB 2 *0 ' (v" cos 0 + w" sin 0) - [EC^" + ECx* (w" cos 0 - v" sin 0) ] ' 60 
+ [EC 1 0” + EC^fw" cos 0 - v" sin 0)]60'loj 


.R 


(64) 


In the equations as well as the boundary conditions, only the linear 
expansion of cos (0 + 0) and sin(0 + 0) should be used to avoid complicating 
the equations with unnecessary small terms. The single underlined terms are 
nonlinear terms and the double underlined terms are the third-order linear 
terms discussed earlier. The tension-torsion term 


30 


- [eak a v(u' * ^ * *£)] 

may appear to be a nonlinear term. However, there is a third-order linear 
approximation based on equations (61a) and (62) such that 

2 2 

EAk^ 2 c()'^u' + + ~ 2 ~ ^ = ^ 2 4>' mft 2 x dx 

+ third- (and higher) order nonlinear terms 


The equations of motion may be solved by Galerkin's method provided that 
natural boundary condition terms in equation (64) that are not identically 
satisfied by the assumed modal functions are added to the resulting modal 
equations in the standard manner (ref. 11) . This operation may be necessary 
to insure symmetric structural and inertial operators. Thus, only the geomet- 
ric boundary conditions (clamped end conditions for a cantilevered beam, e.g.) 
must be satisfied by the assumed modes to obtain equations equivalent to the 
Rayleigh-Ritz method (which deals directly with the energy expressions) . By 
incorporating aerodynamic forces into the equations, stability analyses can be 
performed. 

It is convenient to eliminate u' + (v' 2 /2) + (w ,2 /2) and u from the 6v, 
6w, and 6<|> equations and the boundary conditions before applying a modal 
solution procedure. From the definition of T in equation (62), 
u' + (v' 2 / 2) + (w ' 2 / 2) can be expressed in terms T, v, w, and <t>, while T 
may be determined from equation (61a), thus expressing u' + (v ,2 /2) + (w' 2 /2) 
in terms of v, w, and <p only. Solving for u', integrating over x, and 
then differentiating with respect to t gives an expression for u that may 
be substituted. in equation (61b). The 6u contribution to equation (64) must 
be used in determining u and T; u = 0 at x = 0 and T = 0 at x = R. Thus, 
u may be eliminated from the system of equations. T and T' in equa- 
tions (61b) and (61c) are replaced by equation (61a) for T' and the 
expression 


T = /x R t m ^ 2x 


+ 2£2v) + L u ]dx 


(65) 


for T obtained by integrating equation (61a). Thus, there are now three 
integropartial differential equations in terms of v, w, <J>, and modal func- 
tions are only required for these three deformations. This procedure for 
eliminating u and T is applied to simpler forms of the equations of motion 
in references 4, 8, and 9. 

The equations of motion contain the same basic linear terms as those of 
reference 1 as well as several additional nonlinear terms. A brief discussion 
of the terms should suffice to familiarize the reader with the equations. 
Beginning with equation (65) , the tension is composed of the centrifugal force 
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terms m(fi 2 x + 2ftv) and the applied load L u which is often neglected for 
helicopter applications. The 2mfiv and L u terms are underlined because they 
appear as nonlinear terms in the bending equations (6v and 6w) through the 
centrifugal coupling terms (Tv')' and (Tw')'. The bracketed terms ({ }") in 
the bending equations are bending moments. The 




/ v' 2 w . 2 \ cos ( 0+( » 

V 2 2 /sin(0+4>) 


-EA e A | 

.. Sv 

terms (in the ^ 

equation) 

contain T" as do 


-me 

r«2x cos [Tti + 2£iv°? s 

L sin(0+(})) sin 0J 


and physically they show that the bending moment due to tension is propor- 
tional to e - eA, the distance from the mass centroid axis to the tension 

axis. The -EB 2 *0 '<j> ' C0S ® and +ECi*<f>"^” ® terms represent special twist and 

warp effect terms. The last terms in the brackets are the conventional bend- 
ing moment terms modified in the sense that the cross section is inclined at 
the angle 0 + <J>. The coupling of 4> with v" and w" forms the counterpart 
terms of the product of curvature terms in the torsion equation identified by 
Mil' et at. (ref. 5). To help identify them, these terms are listed sepa- 
rately below. Rearranging the terms with trigonometric substitutions shows 
that the structural coupling terms are proportional to the difference of the 
flap and lead- lag bending stiffness: 

6v: | [El z t - (EI Z « - Elyi ) sin 2 (0 + 4>)]v" + j (EI Z . - EIy«)sin 2(0 + $)w'j 

I f ? 

<5<J>: (EI Z i -EIy»)[v"w" cos 20 +j (w" 2 - v ,,2 )sin 2eJ 

Consider now the remaining terms in the Sv equation. The Coriolis 
term 2mfiu is the counterpart to the -2mflv in the 5u equation. This term 
combines with the nonlinear part of (Tv 1 )' to produce Coriolis (gyroscopic) 
coupling with (Tw')' in the 6w equation. With modal functions, these terms 
form an antisymmetric gyroscopic matrix (ref. 8). The acceleration terms mv 
and -me$ sin 0 are produced by translational and rotational accelerations of 
the cross section, respectively. The offset Coriolis term, 

-2me£l(v' cos 0 + w' sin 0), represents the longitudinal velocity of the cross- 
section mass centroid due to rotation of the cross section during bending. 

This term also combines with -{me[ft 2 x cos(0 + <J>) + 2ftv cos 0]}' to give an 
antisymmetric gyroscopic counterpart and a symmetric stiffness counterpart to 
the -{me[Q 2 x sin(0 + <j>) + 2 ftv sin 0]}' term in the 6w equation. The 
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-mft 2 [v + e cos(0 + <f>)] term is a lateral centrifugal loading due to offset 
from the elastic axis which passes through the center of rotation. Finally, 
-2mQgp C w is a Coriolis term from negative longitudinal velocity due to flap 
bending in the presence of precone. It is antisymmetric with 2mftBp C v of the 
6w equation. The only terms not yet discussed in the <5w equation are the 
accelerations mw and me$ cos 0 which are identical in origin to their 
counterpart terms in the 6v equation and the negative centrifugal force due 
to precone, -mn 2 gp C x. 

The first term in the torsion equation is the tension-torsion coupling 
term, which arises from the tendency of a centrifugal force to untwist a pre- 
twisted blade. The terms involving Bj* and B 2 * are special twist-effect 
terms, which are included in the analysis of reference 1. The terms 

( v 1 2 w' 2 \ 

-EAe^lu' + — + ~ 2 ~J C w " cos 0 - v" sin 0) 

and meft 2 x(w' cos 0 - v' sin 0) are approximately -TeAw" and T'ew' for zero 
0 and thus create a twisting moment due to tension acting on a deflected beam 
at a point offset from the elastic axis. The well-known St. Venant torsion 
rigidity term is -(GJ<j>')'. The [EC^" + EC^w" cos 0 - v" sin 0)]" terms 
arise from the longitudinal stress doing work through the warp displacement 
due to torsion and bending, respectively. The next terms were identified in a 
more simplified analysis by Mil' (ref. 5). These terms have been shown to be 
important for rotor blade stability (ref. 8). The torsion inertia mk m 2 $, and 
the "tennis racquet" terms mn 2 <Hkm 2 - k£j)cos 20 and -mfi 2 (km - k£ )cos 6 sin 6, 
are well-known terms that tend to untwist (or, equivalently, stiffen) the 
rotating blade. The remaining acceleration and centrifugal terms create 
twisting moments because they are offset from the elastic axis. 

These equations in modal form, when linearized with respect to equilib- 
rium for L u = L v = L w = M<j, = 0 , have symmetric stiffness and mass matrices 
and an antisymmetric gyroscopic matrix. This property depends on the reten- 
tion of only those terms from the substitution of the tension T, that are 
consistent with the ordering scheme. 
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PART II 


\ 


DERIVATION OF EQUATIONS BY NEWTON'S SECOND LAW 


In this section, the equations of motion are derived from Newtonian 
mechanics. The derivation treats the forces and moments (elastic, inertial, 
or aerodynamic) applied to the rotating beam. We first consider the equilib- 
rium of a deformed blade in terms of its geometry and the resultant forces and 
moments acting on it. Second, the inertial loadings on the blade are evalu- 
ated in terms of the blade deformation. Finally, the resultant structural 
moments and forces are determined in terms of the blade deformation based on 
the stress-strain and strain-displacement relations. Combining these three 
elements yields a set of integropartial differential equations for the longi- 
tudinal (radial), flapping, lead-lag, and torsion deformations of the elastic 
blade. 


Equations of Structural Equilibrium 

By considering a differential of length dx of the deformed beam, forces 
and moments may be summed to establish the equations of equilibrium (see 
fig. 4). For force equilibrium 


9V , 

* 37 dX - 


V + p dx = 0 


or 


9V 

3x 


+ 



0 


For moment equilibrium 


M 


3M 

3x 


dx - M + 


dxi' 


x 



+ q dx = 0 


or, taking the limit as dx -> 0, 


9M 

8x 


-+ 



x V + q = 0 


( 66 ) 


(67) 


where p, q are "external" forces and moments of inertial or aerodynamic 
origin (aerodynamic forces and moments are taken as given) ; V and M are 
equipollent internal structural forces and moments acting at the elastic axis. 
Also note that moments have been summed about a point at x rather than at 
x * dx.^ 

3 The use of dx rather than distance along the deformed blade might be 
cause for concern. However, to a consistent approximation of neglecting 
squares of slopes compared to one they are the same. 
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The vector equations of equilibrium may be written in several component 
forms. Of the various choices possible, the following seems as convenient as 
any; namely, equations (66) and (67) are expressed in components relative to 
the undeformed x, y, z axes. The moment equations, in scalar component form, 
are 


3M X 

3x 


+ 



v y lx + q x = 0 


3 My 

3x 


+ V x H + q y = 0 


9Mz 

3x 


+ V. 


y 



q z = 0 


(68a) 

(68b) 

(68c) 


This is the same as equation (16) of reference 1, except for a difference in 
sign convention for My and qy. 

Similarly, the force equations are 



op S 

+ 

>d 

X 

ii 

o 

(69a) 


3V„ 

-IT* Py- 0 

(69b) 


OP S 

X tf 3 

+ 

N 

II 

O 

(69c) 

and 

tip 

these may be integrated, noting that V v , V 7 , V Y 

(x = r) , y 

are zero at the rotor 


v x * Px dx 

(70a) 


v y - /x R Py d11 

(70b) 


v x ■ Pz 

(70c) 


More precisely, the upper integration limit should be the deformed length of 
the beam rather than the undeformed length R. However, the deformed length 
is of order R(1 + e xx ) where e xx (the longitudinal strain) may be neglected 
for small rotations - for example, (3w/3x) 2 << 1. 


Anticipating the need for four equilibrium equations for the four 
deformations u, v, w, and <J> , the six equilibrium equations (68) and (69) 
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9M^ 

3x 


9 2 M. 


8x z 

3 2 M z 

ax 2 


X + 


Pz + 


Py ' a5T 


y and V z in the process. 
f68b,c) yields 


* <*z) K * 

q x = 0 

(71a) 

( V * S) * 

8q y n 

l£= 0 

(71b) 

b £) * 

9q z 

ax " 0 

(71c) 


The first force relation, equation (69a), remains 

av Y 


+ p =o 
ax *x 


(71d) 


Either equation (71d) or its integrated form, equation (70a), may be used to 
eliminate V x from the above equations. 

There are now four equilibrium equations and seven unknowns. My, M x , M z , 
V x , v, w, and 4>. The applied loads p x , p y , p z , q x , qy, q z due to inertial 
and/or aerodynamic loadings are functions of v, w, (j>, and a longitudinal 
radial deformation u in the x direction. The inertial relations are 
developed in the next section, and the aerodynamic loadings are assumed to be 
known. In addition, four force and moment deformation relations required to 
express M x , M y , M z , V x in terms of u, v, w, <|> are developed in a later 
section. These latter relations are most easily developed in terms ^ of the 
moment and force components in the deformed body axis system, x', y', z'. 
Hence, the equations of equilibrium will be developed for this axis system. 

The necessary moment transformation between the deformed and undeformed 
axis systems was given by equations (2) and (3); an identical transformation 
holds for the force components V x i, V x , etc. Applying the transformation to 
equations (71) yields 


|m x » ^1 - Myi[v’ cos (0 + <j>) +W' sin(0 + <(>)] 

+M z t [v 1 sin(0 + <f>) - w' cos (0 + <j>) ] J + v 'q y 

+v, |mx' v ' + [My . cos(0 + $ + vV) - M z ' sin(0 + <j> + v'w')]| 

+w ' q z + w' Jm x .w' + ^1 - ^y-^[M y . sin(0 + $) +M Z « cos(0 + 4)] j + q x 


= 0 (72a) 
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(72b) 



cos(6 + $ + v'w») - M z , sin(0 + $> + v'w')] 


)' 


+(V x w')' + p z + q y , = 0 


|m x iw'+ ^1 - sin(0 + $) + M z » cos(0 + $)]J - (V x v» ) ' - p y + q z , =0 (72c) 


where 


V 


x 



dx 


Note that ( ) ' is substituted for 3/ex; this should cause no confusion with 
the primes that denote the deformed body axis system. Equation (72a) may be 
simplified further by first rewriting it as 


M 


x' + jv'(M x .v') . . ^ 

-IK) 


2 

My t cos (0 + $ + v 


) ] + [w'(M x .W)» -^M x , J 

■w.)] + " ~2~ ^ M y » sin(0 + |)J 


-{My, O' cos(0 + <{>) + W* sin (0 + <j>) ] - M z , [v' sin (0 + <j>) - w' cos (0 + <f>)]} ' 


^ y ~ ^ M z' sin(0 + <j> + v'w') + w * ^1 - ^ 


M z , cos(0 + <{>) 


I 


+ q x + y,< fy + W ’ q 2 = 0 <7 3 ) 


Consistent with the ordering scheme, the second and third terms may be 
neglected compared to the first. Also, $ may be replaced by <j>, and the 
squares and products of v f and w ? may be ignored with respect to unity in 
equations (72) and (73) . Expanding derivatives and cancelling terms in equa- 
tion (73) then gives 


M x» " M y f t vM cos ( e + <f>) + w" sin(0 + <[>)] + M z? [v ,f sin (6 + <(>)- w" cos(0 + <J0] 

+ q x + v 'q y + w, q z = ° (74) 


Equations (72b) and (72c) may be simplified further since the product of 
torque M x t and a bending slope v ! or w' may be neglected compared to a 
bending moment M x t or M z i. This is physically evident since a point load on 
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the rotor blade will have a smaller associated moment arm (~ chord or 
thickness! for producing torques than for producing bending moments (~ radius} . 
Sis may be shoL formally b? integrating equations (72a-c), multiplying equa- 
tion (72a) by v' or w', and subtracting the result from the equation (72b) o 
(72c) . This calculation eliminates the terms M x »v' and M x iw’ in equa- 
tions (72b) and (72c), respectively. The additional terms introduced by the 
subtraction will be negligible consistent to second order provided that the 
applied torque loading is small compared to the applied bending moments. 


<bc ^ V 

Therefore, equations (72b) and (72c) may be written 

[My t cos(0 + <)>) - M z > sin(0 + <f>)]" + (V x w')’ + P z + 

[My i sin(0 + 40 + M z » cos (0 + <}>)]" - (V x v 1 ) ' Py + 

The linear approximations for sin(0 + <)>) and cos(0 + <(>) 
used in equations (72) : 


<ly = 0 

(75) 

qz = 0 

(76) 

should 

also be 


cos (0 + 40 - cos 0 - 4> sin 0 
sin(0 + 40 = sin 0 + 4* cos 0 


(77) 


In addition to equations (74) -(76) for moment equilibrium an equation 
for V x . is needed if the tension changes significantly, as a function of the 
rotor blade deformation, from the nominal centrifugal force value due to blade 

rotation 


i aetormatiuii , uum me r 

The equations for V x ., V y ,, and V z , are obtained by transforming 

•a _ i ~ m -ifi -r\ rr 


iUtauuu. * uv iv ' ^ “ 

equations (69a- c) to the deformed axis system, giving 




Vyi[v' cos(0 + <f>) + w' sin(9 + <|>)] 

J rR 

x Px dx 


Vv»v’ + 


^_Xli^[Vy, cos(0 + $ + v’w') - v z . sin(0 + $ + v'w’)] =/ x R p y dx 
V x iw' + - ^-)[Vy* sin(0 + $)+V z i cos(0 + $)]=/ xPz dx 


(78a) 

(78b) 

(78c) 


The required equation for V x i is obtained by substituting Vyi 
equations (78b) and (78c) into equation (78a) . After neglecting 
terms consistent to second order, V x i becomes 


and V z i from 
higher-order 
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V X. - f x \ dx + V 'f x R Py dx + w'^ R p z dx 


(79) 


Equat;i o n s (74)- (76) and (79) are the principal expressions for the equilibrium 
of internal and external forces and moments acting on the blade. For rotor 
blade applications, the second and third terms are negligible to second order- 
thus, V x i = V x . ’ 


Inertial Loading 


The present derivation of the inertia loads closely follows appendix B of 
reference 1 except that nonlinear terms are retained in accordance with the 
ordering scheme. The deformed position of a generic point of the rotor blade 
will be expressed with respect to the rotating undeformed, nonpreconed x„, y, 
Z coordinate system. Note that here (but not elsewhere in this report) , u, 
v > Px> Py> P z > q x , qy, q z are measured with respect to a nonpreconed 

rat er than a preconed coordinate system; the nonpreconed coordinate system 
corresponds to that of reference 1 . The deformations and loadings with 
respect to the two coordinate systems are related by the following transforma- 
tion, if the precone angle g pc is small: 


npc 

= u + e pc w 

o 

= Px + 

6 pcPz 

v pc 

= V 

^ x pc 

= + 

^pc^z 

w pc 

= w 

Pzpc 

= p z - 

6 pcPx 

‘f’pc 

= * 

Q Z 

nz, pc 

= - 

6 pc%c 

x pc 

= X 

Py 

y pc 

■ p y ; 

Qy 

J pc 


The inertia terms are derived using the nonpreconed variables. In the final 
equations for the inertial loadings, however, equation (80) will be used to 
transform the forces, moments, and deformations into preconed coordinate 
system consistent with the remainder of this report. 

As a result of equation (80) , the position of a generic point with 
respect to the x p , y^ Z axes (see eqs. (40)) becomes 


x l - x + u - v' (y x - v) - (w' + B pc ) ( z ! - w - g pc x) - A<f>' 

yi = v + y cos(0 + 4>) -5 sin(0 + <f>) 

z l = w + S pc x + yi sin (6 + <)>)+? cos (6 + <|>) 


(81a) 

(81b) 

(81c) 


Substituting equations (81b) and (81c) into equation (81a) yields 


x i = x + u- A*' - v'[n cos (0 + <f>) - x, sin(0 + $)] -(w' + B pc ) [ n sin(0 + <(.)+? cos(0 + <f>)] 
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The time derivatives are 

x x = u - [v' + (w* + 3p C )<t>] [n cos(6 + 40 - C sin(0 + <f>) ] 

+ (v’<j> - w') [n sin(9 + <J0 + C cos(0 + 40] - X4>' 
x x = ii - [v' + w'4> + (w' + 6 pc )4 + (w' - v ' 4>) 4>] [n cos(0 + 4O - 5 sin(0 + <|>)] 

_X$* + {[v' + (w' + 6p C )i34> + v' 4> +v'$ -w')[n sin(e + 40 + C cos(0 + 40]^ 


y 1 =v - 4> [n sin(0 + 40+5 cos(0 + 40] 

y x = v - i 2 [n cosCe + ♦) - ? sin(0 + 40] - 4>[n sin(0 + 40 + 5 cosC9 + 40] 

= w + <f > [ti cos (0 + 40 - ? sin(0 + 40] 

'ii = w - 4> 2 [n sinC0 + 40 + 5 cosCe + 40] + $[n cos(0 + ♦) - C sin ( e + 40] 


(83b) 


(83c) 


The components of acceleration in the i, j, k coordinate system, including 
the usual contributions due to rotation D, are (ref. 1) 


a x = xj - ft 2 xi - 2Qyi 

ay = yi - fl 2 yr + 2ft*! 

a z = zi 


(84) 


The inertial forces p 1 and moments q 1 may now be evaluated in a 
D’Alembert sense using equations (81) -(84) and the following definitions 


Px 1 = "// A P a x dr > d ? 

Py 1 = - /4 pa y dr) d ^ 

Pz 1 E "// A P a z dn d? 


(85) 


qx 1 5 // A p[ ayCzj - W) - a z ( yi - V)]dn d? 
qy 1 = '// A P a xCzi - «)dn d? 
q z X = // A Pa x (yi - v ) dri d? 


Note that the present sign convention for q * is d *^ e *®" t f g5? m an^higher 1 
Equations (81)- (84) are substituted into equations (85) and (86), and h g 
order terms are discarded. In reference 1, it is assumed that 


40 


cos(0 + <{>) = cos 0 
sin(0 + cj>) = sin 0 


and all products in u, v, w, and their derivatives are ignored - that is, 
equations (85) and (86) are linearized. This is not the same as simply 
assuming <J>, v', etc., are small and using the small-angle approximation. 

Here we apply a consistent set of simplifying approximations to equations (85) 
and (86). In ordering the various terms it is noted that (e/R) 2 , (c/R) 2 , 
(t/R) 2 are small compared to one; hence, products of these quantities with 
themselves or with v/R, w/R or <|> can be neglected with respect to unity. 

The order of magnitude estimates ( )' - 1/R, (•) ~ Q are also used. 

The evaluation of p* and q* from equations (85) and (86) and the 
transformation from nonpreconed to preconed axes leads to the same inertial 
terms that were developed in Part I. To second order, the results for the 
nonpreconed variables are 


Px 1 = m(fl 2 x + 2 Civ) 

Py* = -m{v - e'<p sin 0 - ft 2 [v + e cos(0 + 4>) ] 

+2fi[u - e(v' cos 0 + w' sin 0)]} 

T 

p z = -m(w + e<f> cos 0) 


(87) 


q x J = m(e[(v - ft 2 v)sin(0 + <J>) - w cos(0 + <f>) + 2fiu sin 0] 

a 2 

-k m 2 $ - fi 2 (% 2 - k mi )cos(0 + <J>)sin(0 + <(>) 

- 2fi [( k m 2 " k mi)v' sin 0 cos 0 + w‘( k m 2 sin 2 0 + k^ cos 2 0)]} > 
qy 1 = me[ft 2 x sin(0 + <)>) + 2ftv sin 0] 
q z J = -me[ft 2 x cos (0 + <)>)+ 2ftv cos 0] 


( 88 ) 


Everywhere in equations (87) and (88) , the usual linearized expressions for 
sin(0 + <{>), cos (0 + <f) should be used. Additional simplifications will be 
made when these terms are compared to the elastic forces. For example, the 

third-order terms of q x ^ + v'qy^ + w'q z ^ are neglected except the double 
underlined ones. 


The correspondence between the present inertial forces and those of 
Part I is shown by 

^x 1 + v ’ ( iy I + w ’qz T ) 

Pz 1 + (Qy^ ) 1 < > Z w - Z w . > (89) 

Py 1 - (q 2 X ) ' *■* Zy - Zy. ' ) 
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Force-Deformation and Moment -Deformat ion Relations 

The force and moment resultants are expressed in terms of the 
deformations by resolving the distributed stresses into a resultant force and 
moment system acting at the elastic axis. These resultants are essentially 
the same as those obtained in Part I. Here, they are derived from a physical 
point of view, whereas in the variational method they arise as mathematical 
groups of terms. The resultant axial force is 


V x t = Sf^ xx dn dc = // A Ee xx dn dc 

= eaJu' - e A [v" cos (0 + <f>) + w" sin (9 + 4>)]| (90) 


The bending moments are 


Myi = /J^£cr xx dri d C - «f/ A E£e xx dn d£ 

= Elyi [v" sin(0 + «j>) - w" cos(0 + <j>)] - ECi*^’ 


M z . = -/4na xx dn d? = // A Ens xx dn dc 

2 £ v 

= -EAe A (u' + + EB 2 *6'<}>' + EI Z » [v" cos (0 + <|)) + w" sin(0 + 40] (92) 


The twisting moment with respect to the x' axis is more complicated because 
of the interaction of longitudinal and shear strains with the warp. Brunei le 
(ref. 12) has identified the twisting moment by use of the variational method. 
For the assumptions of the present work, his equation takes on a slightly 
different form 


M 


x> 


•jn 


n° xi - ~ S°xn + ^ 


(^f ' JV°xx ^ ac 


= S X I - (P X |) ' 


(93) 


where 


S x . = // A [(n - A^)o x? - (c + V°xn ]dri = GJr 

p x , e -// A Xa xx dn dc = EC^" + EC 1 *(W” cos 0 - v" sin 0)^ 


(94) 
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The contribution of the P x , term is clearly related to the work done by the 
longitudinal stress acting through the virtual displacement in the longitudi- 
nal direction due to warp. Moreover, as shown in reference 1 (based on a more 
detailed discussion in ref. 13), there is an additional twisting moment due to 
longitudinal stress: 


T X T = (0 + <|>) ' // A (n 2 + s 2 )a xx dn dc 


/ 2 . 2 

= EAk A 2 (0 + J)'(u' + + EB^B'V - EB 2 *9 ' (v" cos e + w' 


sin 0) (95) 


This expression represents the component of force du e to longitudinal stress 
normal to the elastic axis with moment arm A 2 + C 2 . The total resultant 
twisting moment is now expressed by combining equations (93) and (95) : 

M x * = T X ' + S x i - (P x » ) ' = GJ4> ' + EAk A 2 (0 + $) 1 ^u' + ^ 

+EB!*0' V - EB 2 *0'(v" cos 0 + w" sin 0) 

-[ECi<f>" + EC 1 *(w" cos 0 - v" sin 0)]' rg< 


The double underlined terms may be neglected according to the ordering scheme. 
However, as discussed above, these terms are important for configurations with 
low torsion stiffness. The various section constants in the above resultants 
are defined in Part I, equation (35). 

Equations (90) -(96) give the desired force- and moment-deformation 
relations. They can be simplified by making the usual small-angle assumptions: 

cos (0 + <(>)= cos 0 — 4> sin 0 

sin(0 + <)>) = sin 0 + </> cos 6 

In general, no additional approximations can be made, although for special 
blade cross sections, several of the blade section constants may reduce to 
zero - for example, for a doubly symmetric cross section e A = B 2 * = 0. 


Final Equations and Boundary Conditions 

The four equations of equilibrium, the inertial loadings, and the force- 
and moment-deformation relations may be combined to give the final equations 
of motion in terms of the deformations u, v, w, <f>. That is, equations (87), 
(88), transformed to the preconed axes, (90), (91), (92), and (96) are sub- 
stituted into equations (74) -(76) and (79). The resulting equations are sim- 
plified consistent to second order by discarding higher-order terms in 
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accordance with the ordering scheme. The final equations are identical to 
those obtained by the variational method in Part I, and are not repeated here 


As shown in Part I, the boundary conditions for the variational method 
were obtained as a natural by-product of the derivation of the equations of 
motion, without requiring any independent consideration. For the Newtonian 
method, special attention must be given to the boundary conditions, which are 
of two types. The geometric conditions, normally applied at x - 0, are 
simply 

u = v = w = <j> = v'=w’ = 0 for x = 0 (97) 


The natural conditions for the free end x - R may be expressed as 
or, since V x = V x ' to second order, 

M x i = M y . = M z > = V x . = V y = V z = 0 for x = R 


M = V = 0 


(98) 


In addition, the work done by the longitudinal stresses acting through the 
warp displacements must be zero at x = R. Thus 

P„, = 0 for x = R (") 


Restrained warping at the hub implies that 

4 > ’ = o for x = 0 


( 100 ) 


That is, there can be no warp displacement at x = 0. These conditions reduce 
to the standard cantilever boundary conditions when warp effects (C x and Li ) 
are neglected. Recall from equations (68) and (69) that the shears V y , V z ~ 
may be expressed in terms of V x ., M y t, M z ., q y , and q z . When x = R, x ' - 
and the shears may be written (for q x << q y , q z ) 


V y = -M’ sin(0 + <(>) - M z t cos (6 + <j>) - q z = 0 1 
y y “ l for x = R 

V z = M^ , cos (0 + d>) - M z . sin (0 + <f>) + q y = 0 ) 


( 101 ) 


The boundary conditions (97)- (101), when expressed in terms of the deforma 
tions u, v, w, and <J>, are equivalent to equation (64) - the boundary 
conditions as derived from the variational method. 


CONCLUDING REMARKS 


Nonlinear equations of motion for the elastic bending and torsion of 
twisted nonuniform rotor blades have been derived by two complementary methods. 
When deriving a system of nonlinear equations it is advantageous Jo check the 
resulting equations in this manner. The use of the Newtonian method alone, 
for example, may give rise to nonsymmetric structural or inertial operators. 
However, with the variational method, a consistent set of equations will 
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automatically result if the energy expressions are accurate to the desired 
order of magnitude. With the Newtonian method, the boundary conditions must 
be established on the basis of physical reasoning. The variational method 
leads to the appropriate boundary conditions as a by-product of the derivation 
of the equations of motion. The two methods are complementary in* that the 
Newtonian method provides clearer physical understanding of the force compo- 
nents and reactions, while the variational method is more precise mathemati- 
cally - that is, there is less chance of inadvertently leaving out important 
terms or including unnecessary ones. 

In the resulting system of equations, several important nonlinear terms 
are identified. First, the centrifugal term proportional to lead- lag velocity 
in the tension equation combines with the centrifugal coupling terms in the 
bending equations to produce nonlinear flap- lag inertial terms. The longitu- 
dinal velocity in the lead- lag equation, a Coriolis term, is expressed in 
terms of bending quantities as another nonlinear flap-lag inertial term. 

These terms, when linearized with respect to equilibrium, are antisymmetric 
gyroscopic terms and significantly influence hingeless rotor stability and 
forced response phenomena. Second, the nonlinear bending- torsion coupling 
term in the torsion equation is written in a form similar to the one identi- 
fied by Mil' in reference 5. The twisting moment arises from bending in two 
directions and is proportional to the difference in bending stiffness and the 
product of curvatures. The counterpart nonlinear bending- tors ion coupling 
terms in the bending equations appear in the form of a change in elastic 
coupling due to elastic twist. These bending-torsion coupling terms are also 
important in determining the aeroelastic stability of hingeless rotors. 
References 8 and 9 give a preliminary discussion of the application of the 
present equations to a stability analysis of uniform blades. 

The effect of warp interacting with the longitudinal stress is included 
for completeness, although it can be neglected without essential loss of 
accuracy for most applications involving closed cross sections. 

An important step in the derivation of the equation is the development of 
a nonlinear strain-displacement relation based on an exact coordinate trans- 
formation from the undeformed system to the deformed system. This general 
relation between strain and displacement may be useful for extending the equa- 
tions to include higher-order terms than these considered herein. 

Future improvements in the equations should include the incorporation of 
torque offset, root-pitching motion with elastic restraint, droop, pre-lag, 
and other parameters (discussed in ref. 9) to describe more completely the 
behavior of hingeless rotor blades. With the addition of aerodynamic forces, 
the equations of motion may be used for stability analyses and frequency 
response studies of hingeless helicopter rotor blades. 


Ames Research Center 

National Aeronautics and Space Administration 
and 

U. S. Army Air Mobility R§D Laboratory 

Moffett Field, California 94035, Aug. 20, 1974 
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APPENDIX 


DEFORMED BLADE COORDINATE SYSTEM TRANSFORMATION 
D. A. Peters 

In the section devoted to the development of the strain-displacement 
relations, the transformation [T] , given to second order in equations (1), is 
of central importance. Here [T] is determined in terms of u, v, w, and $ 
based upon the development of reference 14 [T] is the transformation between 

the undeformed blade coordinate system 1 , j, k, and the deformed blade system 

i', j', k' : 


(Al) 


\ i 


i’ 

ll 

1, 1 


< — 1 . 

■“V 

II 

j } 

l> 1 


Ik’ I 

kl 

\ / y 

l / 


In terms of Euler angles (fig. 5), [T] may be expressed as 


[T] = 


cos B cos £ 


cos B sin C 


sin B 


-sin 0 sin B cos X, cos 0 cos t, cos B sin 0 

-cos J sin "c -sin ? sin B sin 0 

-cos 0” sin B” cos 5 -sin 0 cos £ cos B cos 0 

sin O' sin T -sin t, sin B cos 6 


(A2) 


The Euler angles are taken in the order B, 0. They uniquely define the 
orientation of the blade principal axes with respect to the undeformed coor- 
dinate system. Rotor bl_ad£ equations are normally written in terms of_ v, w_j_ 
and <|>, rather than B, 9. As shown in figure 5, the Euler angles C and B 
can be easily expressed in terms of v and w, yielding 


[T] = 


A _ v+ 2 - w +2 

v + 

w + 

r — 

- [sin 0w + /l - v +2 - w +2 

[cos 0/1 - v +2 - w +2 

sin 0/1 

+ cos 0v + ]//l - w +2 

-sin 0v + w + ]//l - w +2 


-[cos 0w + /l - v +2 - w +2 

-[sin ?/l - v +2 - w +2 

cos 0/T~ 

+ sin 0v + ]/ /l - w +2 

+cos ?v + w + ]//l - w +2 



• W 


■ w" 


(A3) 
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Figure 5.- Deformation and Euler angles. 

The determination of the third Euler angle 0, however, requires the formula- 
tion and solution of a differential equation for [T] . Consider the small 
rotation wdr of the blade-fixed coordinate system, which occurs as r goes 
through the increment dr. (To second order there is no distinction between 
distance along the deformed dr and undeformed dx elastic axes; hence, 
dr = dx.) The vector components of the rate of rotation u can be identified 
as the torsional rotation rate, = (0 + <j>) + ,l and the bending curvature oi-j 
and uij,. as shown in figure 6. This infinitesimal rotation can be written in 
terms of the transformation matrix [T] , giving the differential equation 

iNote that in this purely geometric exercise there is no distinction 
between pretwist 0 and elastic twist <J>. 
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[T] + = 


0 

-“k 


u“j 


U)R “Wj 

0 

-o)j 0 


[T] 


which can be solved by applying the identity [T] -1 = [T]^: 


'0 

-“k 


.“j 


OJR 

0 

-on 


-W-l 


CDi 


= [T] + [T] T 


(A4) 


Formally, there are nine equations; however, there are only three independent 
ones due to the usual orthonormality conditions on the transformation matrix. 
Substitution of equation (A3) into (A4) , eliminating wj and ojr and identi- 
fying u>i = (9 + <J>) + , yields the exact solution 


e - e + * - r r - dr r 

*n /l - v +2 - w +2 »n 


w+ 2 w ++ v + dr 


(1 - w +2 )/l - v +2 - w +T 


(A5) 


where the pretwist and twist are both taken to be zero at the root. To 
second order, equation (AS) may be written 


9 s Q + (j, _ J V "w* dx 


f q v r» (A6) 

Thus, from (A6) and (A3), [T] may be written to second order as 


[T] = 


1 - 


v- 2 W 2 


W' 


-[v* cos(9+<{>) +w' sin(9+<J>)] cos (9+$+v'w ( ) ^1 - ^ sin(9+$) ^1 - ^ 

2 2 

[v' sin(9 + $) - w' cos (9 + <j>) ] -sin(9+$+v'w’) ^1 - ~ 2 ~) cos (9+$) ^1 - ~ 2 ~ ^ 


(A7) 


where 


<j> = <|> - ,/^ X v"w' dx 


Equation (A7) is the desired result, compare equations (1). 
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It is also useful to record here expressions for coj and which can be 
obtained directly from (A3) and (A4) : 



To second order, these may be simplified as 

o»j = v" sin(0 + <f) - w" cos (6 + <f) 
0 )^ = v" cos (0 + 4>) + w" sin(0 + 4>) 


(A9) 
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